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An analysis of the neutrino oscillation physics capability of a low energy Neutrino Factory is 
presented, including a first simulation of the detector efficiency and event energy threshold. The 
sensitivity of the physics reach to the presence of backgrounds is also studied. We consider a 
representative baseline of 1480 km, we use muons with 4.12 GeV energy and we exploit a very 
conservative estimate of the energy resolution of the detector. Our analysis suggests an impressive 
physics reach for this setup, which can eliminate degenerate solutions, for both large and small 
values of the mixing angle #13, and can determine leptonic CP violation and the neutrino mass 
hierarchy with extraordinary sensitivity. 
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I. INTRODUCTION 

In recent years compelling evidence for neutrino oscillations has been found in experiments with atmospheric [l[ , 
solar 0, [H, EL B S 0]i reactor Q and long-baseline accelerator neutrinos [§, [HI] • Two mass squared differences, 
Am|j = m 2 — m 2 , have been measured with good accuracy, their present best fit values being lAm 2 ^ = 2.5 x 1CP 3 eV 2 

and lAm^xl = 8.0 x 10~ 5 eV 2 . In addition, explaining the experimental data in terms of neutrino oscillations requires 
two large mi xing angles in the leptonic mixing matrix U . Their best fit values are: sin 2 #12 = 0.30 and sin 2 2#23 = 1, 
see Refs. [HI, E2 EH . Despite the remarkable recent progress in our understanding of neutrino physics, fundamental 
questions remain unanswered. It is crucial to establish the nature of neutrinos - whether they are Dirac or Majorana 
particles, the neutrino mass ordering, the absolute neutrino mass scale, the value of the unknown mixing angle #13, 
the presence or absence of CP violation in the leptonic sector, and the precise values of the already known oscillation 
parameters. This information will help shed light on the physics beyond the Standard Model responsible for neutrino 
masses and for the leptonic mixing structure. 

In order to achieve these goals, very sensitive neutrino experiments will be required. In particular, long baseline 
oscillation experiments are expected to play an important role in providing precision measurements of the neutrino 
oscillation parameters, the CP-violating phase S, and a determination of the neutrino mass ordering. 

Neutrino Factories [14j, in which a neutrino beam is generated from muons decaying within the straight sections 
of a storage ring, have been studied extensively in the pas t , a nd have been shown to be sensitive tools for studying 
neutrino oscillation physics [3, ES El EE El, El HI [2l|,j2i, |2l, [11 [11, [11 . In a Neutrino Factory far detector, the 
experimental signature for the so called golden channel |l8| is the presence of a wrong-sign muon [ljj EH; i- e - a m uon 
with opposite sign to the muons stored in the neutrino factory. Wrong-sign muons result from v e — > oscillations, 
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and can be used to measure the unknown mixing angle #13, determine the neutrino mass hierarchy, and search for CP 
violation in the neutrino sector. This physics program requires the detection of charged current (CC) muon-neutrino 
interactions, and the measurement of the sign of the produced muon. If the interacting neutrinos have energies of 
more than a few GeV, standard neutrino detector technology, based on large magnetized sampling calorimeters, can 
be used to measure wrong-sign muons with high efficiency and v ery low backgrounds. This has been shown to work 
for Neutrino Factories with energies of about 20 GeV or greater [3, HH, H] • 

Lower energy Neutrino Factories [III, which store muons with energies < 10 GeV, require a detector technology that 
can detect lower energy muons. Recently ideas have emerged for a Neutrino Factory detector based on a fully active 
calorimeter within a potentially affordable large volume magnet. These ideas encourage consideration of low energy 
Neutrino Factories. Initial studies [2q . based on a first guess for the performance of a low energy Neutrino Factory 
detector, suggested that a Neutrino Factory with an energy of about 4 GeV would enable very precise measurements 
of the neutrino mixing parameters. In the present paper, we consider in more detail the expected performance of a 
low energy Neutrino Factory detector, and update and extend our previous studies to include a more realistic detector 
model and a more comprehensive study of systematic effects. In particular we exploit the low energy threshold of the 
detector and make a very conservative estimate for its energy resolution which, tog ether with the broad spectrum of 
the neutrino factory beam, facilitates the elimination of degeneracies [28l.l29l.l30l.l3ll] It is well known that even a very 
precise measurement of the appearance probability for neutrinos and antineutrinos at a fixed L/E allows different 
solutions for (#13, sign(Am 2 3 ), 6), weakening severely the sensitivity to these parameters. Many strategies have been 
advocated to resolve this issue which in general involve another detector (2fj, |32j, |33j, |3J, |3a, |36|, |37f or the combination 
with another experiment [SIlSllllialiSIMIilSllilliilMI^I- Usin S the energy dependence of the signal 
in the low energy Neutrino Factory, we find that a 4 GeV Neutrino Factory can unambiguously determine all of the 
neutrino oscillation parameters with good precision provided sin 2 2^13 > few -10 -3 . Hence a low energy Neutrino 
Factory would be a precision tool for both large and small 6*13. 

In Section II we describe the design for the low-threshold detector and its performance. In Section III, we discuss 
in detail the physics reach of the proposed setup. We first consider the disappearance signal in order to determine 
precisely the value of the atmospheric mass squared difference and, possibly, the type of hierarchy even for #13 = 0. 
Then, we consider the appearance signals v e — > and v e — > v lll which depend on #13, 5 and the type of neutrino 
mass ordering. We perform a detailed numerical simulation and discuss the sensitivity of the low-energy Neutrino 
Factory to these parameters. In Section IV, we draw our conclusions. 

II. DETECTOR DESIGN AND PERFORMANCE 

A totally active scintillator detector (TASD) has been proposed for a Neutrino Factory, and results from a first 
study of its expected performance are described in the recent International Scoping Study Report [48j] . Using a TASD 
for neutrino physics is not new. Examples are KamLAND Q, which has been operating for several years, and the 
proposed NOi/A detector [Zgj], which is a 15 — 18 Kton liquid scintillator detector that will operate off- axis to the NuMI 
beam line [5(3] at Fermilab. Note that, unlike KamLAND or NO^A, the TASD we are investigating for the low energy 
Neutrino Factory is magnetized and has a segmentation that is approximately 10 times that of NO^A. Magnetization 
of such a large volume (> 30,000 m 3 ) is the main technical challenge in designing a TASD for a Neutrino Factory, 
although R&D to reduce the detector cost (driven in part by the large channel count, 7.5 x 10 6 ) is also needed. 

The Neutrino Factory TASD we are considering consists of long plastic scintillator bars with a triangular cross- 
section arranged in planes which make x and y measurements (we plan to also consider an x-u-v readout scheme). 
Optimization of the cell cross section still needs further study since a true triangular cross section results in tracking 
anomalies at the corners of the triangle. The scintillator bars have a length of 15 m and the triangular cross-section 
has a base of 3 cm and a height of 1.5 cm. We have considered a design using liquid as in NOi/A, but, compared to 
NO^A, the cell size is small (NOi/A uses a 4 x 6 cm 2 cell) and the non-active component due to the PVC extrusions 
that hold the liquid becomes quite large (in NO^A, the scintillator is approximately 70% of the detector mass). Our 
design is an extrapolation of the MINERi/A experiment [5l| which in turn was an extrapolation of the DO preshower 
detectors [52j ]. We are considering a detector mass of approximately 35 Kton (dimensions 15 x 15 x 150 m). We 
believe that an air-core solenoid can produce the field required (0.5 Tesla) to do the physics. 

As was mentioned above, magnetizing the large detector volume presents the main technical challenge for a Neutrino 
Factory TASD. Conventional room temperature magnets are ruled out due to their prohibitive power consumption, 
and conventional superconducting magnets are believed to be too expensive, due to the cost of the enormous cryostats 
needed in a conventional superconducting magnet design. In order to eliminate the cryostat, we have investigated a 
concept based on the superconducting transmission line (STL) that was developed for the Very Large Hadron Collider 
superferric magnets [53[ . The solenoid windings now consist of this superconducting cable which is confined in its own 
cryostat (Fig. [TJ. Each solenoid (10 required for the full detector) consists of 150 turns and requires 7500 m of cable. 
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There is no large vacuum vessel and thus no large vacuum loads which make the cryostats for large conventional 
superconducting magnets very expensive. 

The Neutrino Factory TASD response has been simulated with GEANT4 version 8.1 (Fig.HJ). The GEANT4 model 
of the detector included each of the individual scintillator bars, but did not include edge effects on light collection, or 
the effects of a central wavelength shifting fiber. A uniform 0.5 Tesla magnetic field was simulated. 

Samples of isolated muons in the range of momentum between 100 MeV/c and 15 GeV/c were simulated to allow the 
determination of the momentum resolution and charge identification capabilities. The NUANCE [54j event generator 
was also used to simulate 1 million v e and 1 million interactions. Events were generated in 50 mono-energetic 
neutrino energy bins between 100 MeV and 5 GeV. The results that follow only have one thousand events processed 
through the GEANT4 simulation and reconstruction. 

The detector response was simulated assuming a light yield consistent with MINERVA measurements and current 
photo detector performance [4^ |. In addition, a 2 photo-electron energy resolution was added through Gaussian 
smearing to ensure that the energy resolution used in the following physics analysis would be a worst-case estimate. 
Since a complete pattern recognition algorithm was beyond the scope of our study, for our analysis the Monte 
Carlo information was used to aid in pattern recognition. All digitised hits from a given simulated particle where 
the reconstructed signal was above 0.5 photo electrons were collected. When using the isolated particles, hits in 
neighboring x and y planes were used to determine the 3 dimensional position of the particle. The position resolution 
was found to be approximately 4.5 mm RMS with a central Gaussian of width 2.5 mm . These space points were 
then passed to the RecPack Kalman track fitting package [55l ] . 

For each collection of points, the track fit was performed with an assumed positive and negative charge. The 
momentum resolution and charge misidentification rates were determined by studying the fitted track in each case 
which had the better \ 2 P er degree of freedom. Figure [3] shows the momentum resolution as a function of muon 
momentum. The tracker achieves a resolution of better than 10% over the momentum range studied. Figure 01(a) 
shows the efficiency for reconstructing positive muons as a function of the initial muon momentum. The detector 
becomes fully efficient above 400 MeV. 

The charge mis-identification rate was determined by counting the rate at which the track fit with the incorrect 
charge had a better \ 2 P er degree of freedom than that with the correct charge. Figure 0Jb) shows the charge 
mis-identification rate as a function of the initial muon momentum. 

The neutrino interactions were also reconstructed using the aid of the Monte Carlo information for pattern recogni- 
tion. In an attempt to produce some of the effects of a real pattern recognition algorithm on the detector performance, 
only every fourth hit was collected for track fitting. Tracks were only fit if 10 such hits were found from a given par- 
ticle. The Monte Carlo positions were smeared (Gaussian smearing using the 4.5 mm RMS determined previously) 
and passed to the Kalman track fit. The reconstruction returned: 

• The total momentum vector of all fitted tracks, 

• The momentum vector of the muon (muon ID from MC truth) , 

• The reconstructed and truth energy sum of all the hits that were not in a particle that was fitted, and 

• The reconstructed energy sum of all hits in the event. 

The Vfj, CC event reconstruction efficiency as a function of neutrino energy is shown in Fig. [D(a). The fraction of 
Vfj, CC events with a reconstructed muon is shown in Fig. H^b). In this figure the bands represent the limits of the 
statistical errors for this analysis. 

Based on these initial Neutrino Factory TASD studies, in our phenomenological analysis we assume the detector 
has an effective threshold for measuring muon neutrino CC events at E v = 500 MeV, above which it has an energy 
independent efficiency of 73%. The 73% efficiency is primarily driven by the neutrino interaction kinematics, not by 
the detector tracking efficiency. No charge-ID criterion is applied here. The charge misidentification rate information 
is used as input into the effect of backgrounds on the analysis. 

We note that, to fully understand the backgrounds in the TASD requires a simulation that includes neutrino 
interactions and a full event reconstruction. Although this is beyond the scope of the present study, a consideration of 
backgrounds in the well studied Magnetized Fe-Scintillator detector proposed for the high-energy Neutrino Factory [l8[ 
and The International Scoping Study for a Neutrino Factory [48| motivates the 10~ 3 background (contamination) 
assumption used in this paper for the TASD. Before kinematic cuts, the main backgrounds for the Fe-Scintillator 
detector are muon charge mis-ID, charm decay, pion and kaon decay, and are all of comparable order: 1 — 5 x 10 -4 . 



1 At this stage, the simulation docs not take into account light collection inefficiencies in the corners of the base of the triangle. 
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FIG. 1: Diagram of Superconducting Transmission Line design. 

For the TASD at a low energy Neutrino Factory the muon charge mis-ID rate (Fig. 0Jb)) and the charm decay 
background is suppressed (at the level 4 — 8 x 10 -5 ) due to the low energy beam. Pion and Kaon decay in flight 
become the main background concerns at the 1 — 5 x 10 -4 level. A figure of merit for comparing the TASD to a 
conventional Magnetized Fe-Scintillator detector is the ratio of their respective particle decay lengths to interaction 
lengths. For TASD the ratio is about 1; for the Magnetized Fe-Scintillator detector it is approximately 8. So naively we 
can conclude that the decay background in TASD will be 10 times worse than in the conventional detector ignoring any 
kinematic or topological cuts. However the TASD will have vastly superior kink detection to identify decay- in- flight. 
For example we will typically have 40 hits on the pion track before decay. In addition TASD will have continuous 
dE/dx measurements along the track and better overall energy resolution. We believe that these properties will allow 
us to control backgrounds to the 10 -3 level or better. 

III. PHYSICS REACH OF THE LOW ENERGY NEUTRINO FACTORY 

We have previously mentioned that, by exploiting the energy dependence of the signal, it is possible to extract 
from the measurements the correct values of #13 and S, and eliminate the additional solutions arising from discrete 
ambiguities. In the present study, we include the detector simulation results described in the previous section, which 
suggests a lower energy threshold (500 MeV) than previously assumed (2(|, and an energy resolution dE/E = 30% 2 . 
Above threshold, the detector efficiency for muon neutrino CC events is taken to be 73%. 

In the following we consider the representative baseline L = 1480 km, which corresponds to the distance from 
Fermilab to the Henderson mine. However, we believe that the TASD will not require operation deep underground 
in order to remove backgrounds. Results are similar for other baselines in the 1200-1500 km range. The results are 



2 We have assumed a very conservative dE/E = 30% because at this time the simulation work has not yet produced a number for the 
TASD. Based on NOi/A results, we expect the TASD dE/E to be better than 6% at 2 GeV. 



FIG. 2: Schematic of Totally Active Scintillator Detector. 



presented for the high-statistics scenario described in [26[ as well as for a more aggressive scenario which improves 
the statistics of the old high-statistics scenario by a factor of three, to quantify the benefits of increased detector sizes 
and/or stored-muon luminosities. The high-statistics scenario corresponds to 1 x 10 23 Kton-decays (10 years of data 
taking, with 5 x 10 20 useful muon decays of each sign per year, and a detector fiducial mass times efficiency of 20 Kt). 
The more aggressive scenario corresponds to 3 x 10 23 Kton-decays (which could correspond, for instance, to 10 years 
of data taking, with 1 x 10 21 useful muon decays of each sign per year, and a detector fiducial mass times efficiency 
of 30 Kt). 

Table U shows the number of CC muon events expected in the two scenarios explored here for, respectively, the 
positive and negative muons stored in the Neutrino Factory. Notice that, in the absence of oscillations, there would 
be a few times 10 4 v e CC interactions, which would allow a search for v e — > oscillations with probabilities below 
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statistics 
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13 


22 


25 


11 


(10 23 ) Kt-decays 


3 


39 


66 


77 


34 



TABLE I: Neutrino and antineutrino charged currents interaction rates for L = 14-80 km, for the 10 Kton-decays and the 
3 x 10 23 Kton-decays-statistics scenarios. 

All numerical results reported in the next subsections have been obtained with the exact formulae for the oscillation 
probabilities. Unless specified otherwise, we take the following central values for the remaining oscillation parameters: 
sin 2 0i2 = 0.29, Am 2 ! = 8 x 10~ 5 eV 2 , lAm 2 ^ = 2.5 x 10~ 3 eV 2 and 6 23 = 40°. We show in Tables [H and [ml for two 
representative values of 6*13 = 1° and 8°, and the CP phase S — 0°, 90°, 180° and 270°, the number of wrong-sign muon 
events in the two scenarios explored here, for, respectively, the positive and negative muons stored in the Neutrino 
Factory, for normal (inverted) hierarchy. 

For our analysis, we use the following \ 2 definition 

X 2 = E E - N hP )C~ p \ hp ,(n iy - N jy ) , (3.1) 

i,j P,p' 

where iVj ± is the predicted number of muons for a certain oscillation hypothesis, n, lP are the simulated "data" from 
a Gaussian or Poisson smearing and C is the 2Na n x 2Nun covariance matrix given by: 

C i,p\j, P ' = SijSpp'iSni^) 2 (3.2) 
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FIG. 3: Momentum resolution as a function of the muon momentum. 



statistics (Kt-decays) 


5(°) 


/i + stored (wrong-sign: ^ ) 


(i stored (wrong-sign: /i + ) 


1 x 10 23 





880 (340) 


180 (520) 


90 


1230 (505) 


90 (330) 


180 


1000 (340) 


170 (440) 


270 


645 (175) 


260 (625) 


3 x 10 23 





2640 (1020) 


540 (1550) 


90 


3700 (1520) 


270 (990) 


180 


2990 (1020) 


510 (1310) 


270 


1930 (520) 


780 (1870) 



TABLE II: Wrong sign muon event rates for normal (inverted) hierarchy, assuming v e — + (v e — » v^) oscillations in a 20 Kt 
fiducial volume detector, for a L = 1480 km baseline. We assume here #13 = 8°, i.e. sin 2 28iz ~ 0.076. We present the results 
for several possible values of the CP-violating phase S for both the two scenarios. 

where (Sni :P ) = \Jni^ v + (f sys • n i,p) 2 contains both statistical and a 2% overall systematic error (f sys = 0.02). 

A. Exploring the disappearance channel 

Consider first the disap pea rance channels, already considered in the context of Neutrino Factories [It], and 
carefully explored in Ref. [57|. In Ref. it was shown that, with its high statistics and good energy resolution, 
a low energy neutrino factory can be used to precisely determine the atmospheric neutrino oscillation parameters, 
823 and Amjjp In particular, for an exposure of 3 x 10 22 Kton-decays for each muon sign, and allowing for a 2% 
systematic uncertainty, it was shown that: (i) Maximal mixing in the 23-sector could be excluded at 99% CL if 
sin 2 623 < 0.48 (023 < 43.8°), independently of the value of #13, and (ii) For a large value of #13, i.e. #13 > 8°, the 
#23~octant degeneracy would be resolved at the 99% CL for sin 2 9 2 3 < 0.44 (#23 < 41.5°). In our present study, the 
good energy resolution of the TASD provides sensitivity to the oscillatory pattern of the disappearance signal that is 
comparable to, and somewhat better than, we previously assumed. 
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FIG. 4: (a) Efficiency for reconstructing positive muons. (b) Muon charge mis-identification rate as a function of the initial 
muon momentum. 
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FIG. 5: (a) Reconstruction efficiency of NuMu CC events as a function of neutrino interaction energy, (b) Fraction of NuMu 
CC events with a reconstructed muon. 



In Fig. [S] we show the 68%, 90% and 95% CL contours (for 2 d.o.f) resulting from the fits to the measured energy 
dependent and CC rates at L = 1480 km. Results correspond to 1 x 10 23 Kton-decays, and are shown for Amf^ = 
2.5 x 10~ 3 eV 2 and two simulated values of sin 2 6> 23 (= 0.4 and 0.44). For 13 = 0, P v ^ v „ ($23) = P Vlt -+u„ (tt/2 - #23), 
i.e. the disappearance channel is symmetric under 623 — » tt/2 — 623. However, when a rather large non- vanishing value 
of #13 is switched on, a 023 asymmetry appears in the P v _>j, . Notice that the asymmetry grows with increasing 
#13 and the four- fold degeneracy in the atmospheric neutrino parameters is resolved more easily. We conclude that, 
using only the ^-disappearance data, the uncertainty on Amgi could be reduced down to the 1% — 2% level. In 
principle, the v e disappearance channel could also be used, which is sensitive to #13 and matter effects. However, 
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160 (150) 


80 (110) 
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300 (210) 
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180 


200 (150) 


230 (250) 


270 


65 (90) 


110 (150) 



TABLE III: As Tabled but for 9 13 = 1°, i.e. sin 2 26>i 3 ~ 0.001. 



charge discrimination for electrons has not yet been adequately studied to determine the relevant TASD performance 
parameters. 

The extremely good determination of the atmospheric mass squared difference opens the possibility to determine 
the mass hierarchy by exploiting the effects of the solar mass squared difference on the v, L disappearance probability, 
even for negligible values of 6*13. This strategy was studied in detail in Refs. [58|, H^, [l5(|. The vacuum — > 
oscillation probability is given by 

P(v, -„„) = !- 4|CV| 2 |[/, 2 | 2 sin 2 - A\U^\U^ sin 2 - 4|C/ Al2 | 2 |[/ Al3 | 2 sin 2 ^± , (3.3) 

where the usual notation is used for the mass squared differences Am| and for the elements of the leptonic mixing 
matrix U. In the following we take 6q 3 = 0. The oscillation probabilities depend on whether |Amf 3 | > |Ato2 3 | 
(normal hierarchy) or |Ato 2 3 | < lAm 2 ^ (inverted hierarchy). Precisely measured disappearance probabilities can 
distinguish between normal and inverted hierarchies if there is sensitivity to effects driven by both |Am 2 3 | and Am 2 2 . 
This requires the atmospheric mass squared difference to be measured at different L/E with a precision of better 
than lAm^ |/| Am^ ~ 0.026. In fact, it was pointed out in Ref. [58[ that, for a fixed L/E, the disappearance 
probabilities for the normal and inverted hierarchies are the same if |Ato 2 3 | is substituted with — |Am 2 3 + Am 2 2 + 

^ arctan ^ cos 26>i 2 tan A ™£ L ^ ■ In order to break this degeneracy it is necessary to measure the atmospheric mass 

squared difference at different energies and at distances for which the oscillations driven by the solar term are non 
negligible. In our setup, if we assume a 0% (2%) overall systematic error, we find that the hierarchy can be measured at 
the la level (la level) for the 10 23 Kton-decays case, while for the 3 x 10 23 Kton-decays scenario it can be determined 
at 4cr level (2a level). Note that the systematic errors play a crucial role. It is in principle possible to reduce the 
impact of the systematics errors using the ratios of the number of events at the near and far detectors: 

-m - Ig , (3.4) 

where N-^^p^(i/^[v e \) refer to the number of [D e ] events in the near (far) detector for a fixed energy E. Very good 
energy resolution is required for such cancellations to be effective. In this case, a low energy Neutrino Factory can 
give important information on the type of hierarchy even if 613 = 0. 

B. Simultaneous fits to 6*13 and S 

Next, we study the extraction of the unknown parameters (9i 3 and 8, using the golden channel (u e (p e ) — > v^ip^)). 
We start by considering a neutrino factory scenario with 1 x 10 23 Kton-decays. We find that, for values of 613 > 2°, 
the sign degeneracy is resolved at the 95% C.L. Note that for 613 > 4° the octant degeneracy has already been resolved 
using the disappearance data. 

Figure [7] shows, for a fit to the simulated data at a baseline L = 1480 km, the 68%, 90% and 95% CL contours 
in the (9i3, (5)-plane. Results are shown for background levels set to zero (left panel) and 10~ 3 (right panel) for 
the 10 23 Kton-decays scenario. The four sets of contours correspond to four simulated test points in the (#i 3 ,<5)- 
planc, which are depicted by a star. The simulations are for the normal mass hierarchy and 623 in the first octant 
(sin 2 623 — 0.41 which corresponds to 623 = 40°). Our analysis includes the study of the discrete degeneracies. That 
is, we have fitted the data assuming both the right and wrong hierarchies, and the right and wrong choices for the 
#23 octant. If present, the additional solutions associated to the #23 octant ambiguity are shown as dotted contours. 
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FIG. 6: 68%, 90% and 95% (2 d.o.f) CL contours resulting from the fits at L — 1480 km assuming two central values for 
sin 2 823 = 0.4 and 0.44 and Am§i = 2.5 x 10~ 3 eV 2 . In the left (right) panel, #13 = 8° (3° ). The statistics considered for both 
simulations corresponds to 1 x 10 23 Kton- decays. Only disappearance data have been used to perform these plots. 



Notice from Fig. [7] that the sign ambiguity is resolved at the 95% CL in the 10 23 Kton-decays scenario. Additional 
solutions associated to the wrong choice of the #23 octant are still present in the 10 23 Kton-decays scenario, but notice 
that the presence of these additional solutions does not interfere with a measurement of the CP violating phase 8 and 
013, since the locations of the fake solutions in the ($13, 8) plane, are almost the same as the correct locations. 

The effect of the background can be easily understood in terms of the statistics presented in Tables |TT] and IIIII For 
small values of #13, the addition of the background has a larger impact for 8 ~ —90°, since for that value of the CP 
phase the statistics are dominated by the antineutrino channel, which suffers from a larger background (from z/^'s) 
than the neutrino channel (from C"^'s). For a background level smaller than ~ 10~ 4 , the results are indistinguishable 
from the zero background case. 

We illustrate the corresponding results for the improved scenario of 3 x 10 23 Kton-decays in Fig. [8] Note that the 
higher statistics allow us to consider a smaller value for #13 = 1°. The additional solutions arising from the wrong 
choice for the neutrino mass hierarchy or 6*23 octant are not present at the 95% CL. Furthermore, the addition of a 
background level of 10 -3 does not significantly affect the resolution of the degeneracies, and has only an impact on 
the CP violation measurement. 

The performance of the low energy neutrino factory in the two high statistics scenarios explored here is unique. 
The sign(Am§ 1 ) can be determined at the 95% CL in the 10 23 Kton-decays (3 x 10 23 Kton-decays) scenario if 6*13 > 2° 
(> 1°) for all values of the CP phase <5. The #23-octant ambiguity can be removed at the 95% CL down to roughly 
9 13 > 0.5 — 1.0° for the representative choice of sin 2 623 — 0.41, independently of the value of 8, except for some 
intermediate values of 813 ~ 2° , for which the 623 degeneracy is still present for some values of the CP violating phase 
8. Resolving the #23-octant degeneracy therefore is easier for small values of 613 < 2°. This is due to the fact that, as 
explored in Ref. [261] . the #23-octant degeneracy is resolved using the information from the low energy bins, which are 
sensitive to the solar term. For the setup described in this paper, the solar term starts to be important if #13 < 2°. 
However, notice that the presence of the 623 octant ambiguity at #13 ~ 2° will not interfere with the extraction of 
#13 and 8, since the locations of the degenerate (fake) solutions almost coincide with the positions of "true" , nature 
solutions. 

In Figs. O and [TO] we summarize, for the 10 23 and 3 x 10 23 Kton decays scenarios, the physics reach for a TASD 
detector located 1480 km from a low energy neutrino factory. The analysis takes into account the impact of both the 
intrinsic and discrete degeneracies. Figure [9] shows the region in the (sin 2 2^13, "fraction of 8") plane for which the 
mass hierarchy can be resolved at the 95% CL (1 d.o.f). Contours are shown for zero background, and for when a 
background level of 10 -3 is included in the analysis. Note that, with a background level of ~ 10 -3 , the hierarchy can 
be determined in both scenarios if sin 2 2#i3 > few 10 -3 (i.e. 813 > 2 — 3°) for all values of the CP violating phase 8. 
For a background level smaller than ~ 10 -4 , the results are indistinguishable from the zero background case. 

Figure [TO] shows the region in the (sin 2 26*i 3 , 8) plane for which a given (non-zero) CP violating value of the CP- 
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FIG. 7: 68%, 90% and 95% (2 d.o.f) CL contours resulting from the fits at L — 1480 km assuming four central values for 
S = 0°, 90°, —90° and 180° and #13 = 2° without backgrounds (left panel) and with a background level of 10 -3 (right panel). 
The additional 823 octant solutions are depicted in dotted blue. The statistics considered for both simulations corresponds to 
10 23 Kton- decays. 



phase S can be distinguished at the 95% CL (1 d.o.f) from the CP conserving case, i.e. S = 0, ±180°. The results 
are given for the two statistics scenarios studied here. Note that, even in the presence of a 10 -3 background level, 
the CP violating phase 5 could be measured with a 95% CL precision of better than 20° in the 10 23 Kton-decays 
(3 x 10 23 kton-deacys) luminosity scenario if sin 26»i3 > 0.01 (sin 2 26 13 > 0.002). 

IV. SUMMARY AND CONCLUSIONS 

We have studied the physics reach of a low energy neutrino factory, first presented in Ref. [26] , in which the stored 
muons have an energy of 4.12 GeV. The simulated detector performance is based upon a magnetized Totally Active 
Scintillator Detector. Our simulations suggest this detector will have a threshold for measuring muon neutrino CC 
interactions of about 500 MeV and an energy independent efficiency of about 73% above threshold. We have assumed 
a conservative energy resolution of 30% for the detector. 

In our analysis, we consider the representative baseline of 1480 Km, divide the simulated observed neutrino event 
spectrum into 9 energy bins above the 500 MeV threshold, and exploit both the disappearance (V M — > v^) and the 
golden (y e — > v^) channels by measuring CC events tagged by "right-sign" and "wrong-sign" muons. The results 
can be easily generalized to other baselines in the 1200-1500 km range. We have investigated the dependence of the 
physics sensitivity on statistics by considering a high statistics scenario corresponding to 1 x 10 23 Kton-decays for 
each muon sign, and a more aggressive scenario corresponding to 3 x 10 23 Kton-decays for each muon sign. We have 
also explored the impact of backgrounds to the wrong-sign muon signal by considering background levels of zero and 
10~ 3 . 

We find that, based only on the disappearance channel, maximal atmospheric neutrino mixing can be excluded at 
95% CL if sin 2 6» 23 < 0.44 (0 23 < 41.5°). The atmosphe ric mass difference could be measured with a precision of 
1% — 2%, opening the possibility of determining the neutrino mass hierarchy even if 6*13 = 0, provided systematic 
uncertainties can be controlled. Neglecting systematic uncertainties, the mass hierarchy could be determined at the 
la level (4a level) in the 1 x 10 23 Kton-decays (3 x 10 23 Kton-decays) statistics scenario. 

The rich oscillation pattern of the v e — > (p e — > D^) appearance channels at energies between 0.5 and 4 GeV for 
baselines 0(1000) km facilitates an elimination of the degenerate solutions. If the atmospheric mixing angle is not 
maximal, for the representative choice of sin 2 $23 = 0.4, the octant in which 6*23 lies could be extracted at the 95% 
CL in both scenarios if 6*13 > 0.5 — 1°, for all values of the CP violating phase 5, except for some intermediate values 
of the mixing angle 6*13 in which the fake solutions's location coincides with the true's solution position and therefore 
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FIG. 8: 68%, 90% and 95% (2 d.o.f) CL contours resulting from the fits at L — 1480 km assuming four central values for 
S = 0°, 90°, —90° and 180° and 8n = 1° without backgrounds (left panel) and with a background level of 10 -3 (right panel). 
The statistics considered for both simulations corresponds to 3 x 10 23 Kton-decays. 

the presence of these fake solutions does not interfere with the extraction of 8c p and 6*13. 

In the 10 23 kton-decays scenario, if the background level is ~ 1CP 3 (10~ 4 ), the neutrino mass hierarchy could be 
determined at the 95% CL, and the CP violating phase S could be measured with a 95% CL precision of better than 
20°, if sin 2 29^ > 0.01 (~ 0.006). With a factor of three improvement in the former statistics, the numbers quoted 
above are sin 2 29\3 = 0.005 and sin 2 26*13 = 0.002, for background levels of 10 -3 and 10~ 4 , respectively. In our analysis 
we have included a 2% systematic error on all measured event rates. 

In summary, the low statistics low energy Neutrino Factory scenario we have described, with a background level 
of 10 -3 , for both large and very small values of 6*13 would be able to eliminate ambiguous solutions, determine 613, 
the mass hierarchy, and search for CP violation. Higher statistics and lower backgrounds would further improve the 
sensitivity, and may enable the mass hierarchy to be determined even if #13 = 0. 
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Appendix 

All detector concepts for the Neutrino Factory (NF) require a magnetic field in order to determine the sign of muon 
(or possibly the electron) produced in the neutrino interaction. For the baseline detector, this is done with magnetized 
iron. Technically this is very straightforward, although for the 100 Kton baseline detector does present challenges 
because of its size. The cost of this magnetic solution is felt to be manageable. Magnetic solutions for the TASD 
become much more problematic. The solution that we propose is to use the Superconducting Transmission Line 
developed for the VLHC [53[, Fig.Q] as windings for very large solenoids that form a magnetic cavern (see Fig. [11]) for 
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FIG. 9: 95% CL (1 d.o.f) hierarchy resolution assuming that the far detector is located at a distance o/1480 km at the Henderson 
mine. The solid (dotted) red curves depict the results assuming 1 x 10 23 Kton-decays (3 x 10 23 Kton-decays) without backgrounds. 
The long-dashed (short-dashed) black curves depict the results assuming 1 x 10 23 Kton-decays (3 x 10 23 Kton-decays) with a 
background level of 10~ 3 . 



the detector. The Superconducting Transmission Line (STL) consists of a superconducting cable inside a cryo-pipe 
cooled by supercritical liquid helium at 4.5 — 6.0 K placed inside a co-axial cryostat. It consists of a perforated 
Invar tube, a copper stabilized superconducting cable, an Invar helium pipe, the cold pipe support system, a thermal 
shield covered by multilayer supcrinsulation, and the vacuum shell. One of the possible STL designs developed for 
the VLHC is shown in Fig. [IT] within the main text. The STL is designed to carry a current of 100 kA at 6.5 K in a 
magnetic field up to 1 T. This provides a 50% current margin with respect to the required current in order to reach a 
field of 0.5 T. This operating margin can compensate for temperature variations, mechanical or other perturbations 
in the system. 

The solenoid windings now consist of this superconducting cable which is confined in its own cryostat. Each solenoid 
consists of 150 turns and requires ~ 7500 m of cable. There is no large vacuum vessel and access to the detectors 
can be made through the winding support cylinder since the STL does not need to be close-packed in order to reach 
an acceptable field. We have performed a simulation of the Magnetic Cavern concept using STL solenoids and the 
results are shown in Fig. |T2j With the iron end- walls (1 m thick), the average field in the XZ plane is approximately 
0.58 T at an excitation current of 50 kA. This figure shows the field uniformity in the X-Z plane which is better than 
±2% throughout the majority of the volume with approximately 20% variations near the end-irons. Figure [TBI shows 
the on-axis B field as a function of position along the z axis (in meters). 

We have not yet been able to do a detailed costing of the magnetic cavern. The STL costs can be estimated quite 
accurately (30%) from the VLHC work and current SC cable costs and are believed to be S50M. The total magnetic 
cavern cost is estimated to be less than $150M. This is to be compared to a fully loaded cost savings of the Low-Energy 
Neutrino Factory (compared to the 50 GeV design) as indicated in Ref . [48j . 



[1] Y. Ashie et al. [Super-Kamiokande Collaboration], Phys. Rev. D 71, 112005 (2005). 

[2] B. T. Cleveland et al, Astrophys. J. 496, 505 (1998); Y. Fukuda et al. [Kamiokande Collaboration], Phys. Rev. Lett. 
77, 1683 (1996); J. N. Abdurashitov et al. [SAGE Collaboration], J. Exp. Theor. Phys. 95, 181 (2002); W. Hampel et al. 
[GALLEX Collaboration], Phys. Lett. B 447, 127 (1999); T. A. Kirsten [GNO Collaboration], Nucl. Phys. Proc. Suppl. 
118, 33 (2003). 



13 




FIG. 10: 95% CL (1 d.o.f) CP Violation extraction assuming that the far detector is located at a distance of 1480. The solid 
(dotted) red curves depict the results assuming 1 x 10 23 Kton-decays (3 x 10 23 Kton-decays) without backgrounds. The long- 
dashed (short-dashed) black curves depict the results assuming 1 x 10 23 Kton-decays (3 x 10 23 Kton-decays) with a background 
level o/10" 3 . 



[3] S. Fukuda et al. [Super-Kamiokande Collaboration], Phys. Lett. B 539, 179 (2002). 
[4] Q. R. Ahmad et al. [SNO Collaboration], Phys. Rev. Lett. 87, 071301 (2001). 

[5] Q. R. Ahmad et al. [SNO Collaboration], Phys. Rev. Lett. 89, 011301 (2002) and ibid. 89, 011302 (2002). 

[6] S. N. Ahmed et al. [SNO Collaboration], Phys. Rev. Lett. 92, 181301 (2004). 

[7] B. Aharmim et al. [SNO Collaboration], Phys. Rev. C 72, 055502 (2005). 

[8] K. Eguchi et al. [KamLAND Collaboration], Phys. Rev. Lett. 90, 021802 (2003). 

[9] M. H. Ahn et al. [K2K Collaboration], Phys. Rev. D 74, 072003 (2006). 
[10] D. G. Michael et al. [MINOS Collaboration], Phys. Rev. Lett. 97, 191801 (2006). 
[11] T. Schwetz, Phys. Scripta T127, 1 (2006). 

[12] G. L. Fogli et al, Phys. Rev. D 75, 053001 (2007). 

[13] M. C. Gonzalez-Garcia and M. Maltoni. larXiv:0704.1800l [hep-ph], 

[14] S. Geer, Phys. Rev. D 57, 6989 (1998) [Erratum-ibid. D 59, 039903 (1999)]. 

[15] A. De Rujula, M. B. Gavela and P. Hernandez, Nucl. Phys. B 547, 21 (1999). 

[16] V. D. Barger, S. Geer and K. Whisnant, Phys. Rev. D 61, 053004 (2000); A. Donini et al, Nucl. Phys. B 574, 23 (2000); 

V. D. Barger et al, Phys. Rev. D 62, 073002 (2000). 
[17] V. D. Barger et al, Phys. Rev. D 62, 013004 (2000). 

[18] A. Cervera, et al Nucl. Phys. B 579, 17 (2000) [Erratum-ibid. B 593, 731 (2001)]; 

[19] M. Freund, P. Huber and M. Lindner, Nucl. Phys. B 585, 105 (2000); V. D. Barger et al, Phys. Lett. B 485, 379 (2000); 
J. Burguet-Castell et al, Nucl. Phys. B 608, 301 (2001); M. Freund, P. Huber and M. Lindner, Nucl. Phys. B 615, 331 
(2001). 

[20] A. Donini, D. Meloni and P. Migliozzi, Nucl. Phys. B 646, 321 (2002); 

D. Autiero et al, Eur. Phys. J. C 33, 243 (2004). 
[21] C. Albright et al. , |arXiv:hep-ex/0008064l 
[22] A. Blondel et al, N ucl. Instrum. Meth. A 4 51, 102 (2000); 

M. Apollonio et al, |arXiv:hep-ph/0210192| 
[23] C. Albright et al. [Neutrino Factory/Muon Collider Collaboration], |arXiv:physics/041 1123 
[24] O. Mena, Mod. Phys. Lett. A 20, 1 (2005). 

[25] P. Huber, M. Lindner, M. Rolinec and W. Winter, Phys. Rev. D 74, 073003 (2006). 
[26] S. Geer, O. Mena and S. Pascoli, Phys. Rev . D 75, 093001 (2007). 
[27] P. Huber and W. Winter. larXiv:0706.2862l [hep-ph] . 
[28] G. L. Fogli and E. Lisi, Phys. Rev. D 54, 3667 (1996). 



14 




FIG. 11: Simulation results for magnetic cavern design. 



[29] H. Minakata and H. Nunokawa, JHEP 0110, 001 (2001). 

[30] V. D. Barger, S. Geer, R. Raja and K. Whisnant, Phys. Rev. D 63, 113011 (2001). 

[31] T. Kajita, H. Minakata and H. Nunokawa, Phys. Lett. B 528, 245 (2002); H. Minakata, H. Nunokawa and S. J. Parke, 
Phys. Rev. D 66, 093012 (2002); P. Huber, M. Lindner and W. Winter, Nucl. Phys. B 645, 3 (2002); A. Donini, D. Meloni 
and S. Rigolin, JHEP 0406, 011 (2004); M. Aoki, K. Hagiwara and N. Okamura, Phys. Lett. B 606, 371 (2005); O. Yasuda, 
New J. Phys. 6, 83 (2004); O. Mena and S. J. Parke, Phys. Rev. D 72, 053003 (2005). 
H. Minakata and H. Nunokawa, Phys. Lett. B 413, 369 (1997). 
V. Barger, D. Marfatia and K. Whisnant, Phys. Rev. D 66, 053007 (2002). 
O. Mena Requejo, S. Palomares-Ruiz and S. Pascoli, Phys. Rev. D 72, 053002 (2005). 

M. Ishitsuka et al, Phys. Rev. D 72, 033003 (2005); K. Hagiwara, N. Okamura and K. i. Senda, Phys. Lett. B 637, 266 
(2006). 

O. Mena, S. Pa lomares-Ruiz and S. Pa scoli, Phys. Rev. D 73, 073007 (2006). 
T. Kajita et al, arXiv:hep-ph/ 0609286 
J. Burguet-Castell et al, Nucl. Phys. B 646, 301 (2002); 
P. Huber, M. Lindner and W. Winter, Nucl. Phys. B 654, 3 (2003). 
H. Minakata, H. Nunokawa and S. J. Parke, Phys. Rev. D 68, 013010 (2003). 
V. Barger, D. Marfatia and K. Whisnant, Phys. Lett. B 560, 75 (2003). 

K. Whisnant, J. M. Yang and B. L. Young, Phys. Rev. D 67, 013004 (2003); P. Huber et al, Nucl. Phys. B 665, 487 
(2003); P. Huber et al, Phys. Rev. D 70, 073014 (2004); A. Donini, E. Fernandez-Martinez and S. Rigolin, Phys. Lett. B 
621, 276 (2005). 

O. Mena and S. J. Parke, Phys. Rev. D 70, 093011 (2004). 
P. Huber, M. Maltoni and T. Schwetz, Phys. Rev. D 71, 053006 (2005). 
S. Choubey and P. Roy, Phys. Rev. D 73, 013006 (2006). 

A. Blondel et al, Acta Phys. Polon. B 37, 2077 (2006). A. Blondel lhep-ph/0606 1 lT] 

O. Mena, H. Nunokawa and S. J. Parke, Phys. Rev. D 75, 033002 (2007). O. Mena, |hep-ph/0609031| 
The Internatio nal Scoping St udy for a Neutrino Factory RAL-TR-2007-24 
D. Ayres et al, hep-ex/0503053 
J. Hylen et al, FERMILAB-TM-2018, (1997). 
K. S. McFarland, Eur. Phys. J. A 24S2, 187 (2005). 



15 



l til >■:: .■■>-* ,h 



7 i<:zi-c-an 




\F" VECTOR FIELDS 



FIG. 12: STL Solenoid Magnetic Cavern Field Uniformity in XZ plane. 
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FIG. 13: On-axis B field in T as a function of position along the z axis (m). 
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